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Abstract 

We investigate quasi-Monte Carlo rules for the numerical integration of multivari¬ 
ate periodic functions from Besov spaces S/ q B(T d ) with dominating mixed smoothness 
1/p < r < 2. We show that order 2 digital nets achieve the optimal rate of convergence 
N~ r (\og jVjU- 1 )! 1-1 /?). The logarithmic term does not depend on r and hence improves 
the known bound of Dick [B] for the special case of Sobolev spaces H/ nix (T d ). Secondly, the 
rate of convergence is independent of the integrability p of the Besov space, which allows 
for sacrificing integrability in order to gain Besov regularity. Our method combines char¬ 
acterizations of periodic Besov spaces with dominating mixed smoothness via Faber bases 
with sharp estimates of Haar coefficients for the discrepancy function of order 2 digital 
nets. Moreover, we provide numerical computations which indicate that this bound also 
holds for the case r = 2. 


1 Introduction 

Quasi-Monte Carlo methods play an important role for the efficient numerical integration of 
multivariate functions. Many real world problems, for instance, from finance, quantum physics, 
meteorology, etc., require the computation of integrals of d -variate functions where d may be 
very large. This can almost never be done analytically since often the available information 
of the signal or function / is highly incomplete or simply no closed-form solution exists. A 
quasi-Monte Carlo rule approximates the integral 1(f) = J^ 01 ^ d f(x)dx by (deterministically) 
averaging over N function values taken at fixed points X,v = {a; 1 ,..., x N }, i.e., 

1 N 

lN(X N J):=-J2f (**), 

V i =1 

where the d-variate function / is assumed to belong to some (quasi-)normed function space 
Fd C C([0, l] rf ). Since the integration weights are positive and sum up to 1, QMC integration 
is stable and easy to implement which significantly contributed to its popularity. The QMC- 
optimal worst-case error with respect to the class Fd is given by 

QMC N (F d ):= inf sup \I(f) - I N (X N , f)\. (f.l) 

*ArC[0,l] d ||/|F d ||<l 
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In this paper we investigate the asymptotical properties of Dick’s construction [6| of order a 
digital nets Xjy where N = 2 n . This construction has recently attracted much attention in 
the area of uncertainty quantification mi mg. In the present paper we are interested in the 


asymptotic optimality of those higher order nets in the sense of (1.1) with respect to being 
a periodic Nikol’skij-Besov space S/ q B(T d ) with smoothness r larger than 1 and less than 2. 
Dick [6j showed for periodic Sobolev spaces i x (T d ) = S , 2 2 i?(T fZ ^ 


QMC N (H r mix (T d )) < N-^(\ogN) d ^ 


rj — 1 


N > 2. 


( 1 . 2 ) 


if 1/2 < r < a. He also considered non-periodic integrands, see [7]. However, well-known 
asymptotically optimal results for the integration of periodic Sobolev functions, see for instance 
the survey [50] . show that the exponent of the log should be independent of the smoothness 
parameter r, namely (d — l)/2. In that sense, 0 is far from being optimal. Nevertheless, 
Dick’s bound (1.2) beats the well-known sparse grid bound if r is an integer and d is large. 


see 


and (1.8) below, which 


The latter bound involves the log-term (log jV)( d_1 '(r+i/ 2 )^ 
represents the best possible rate among all cubature formulas taking function values on a sparse 
grid [20] . 

The aim of this paper is twofold. On the one hand we aim at showing the sharp relation 


Q,MC N (H/ aix (T d )) x N~ r (log N)^ 1 ^ 2 , N > 2 


(1.3) 


if 1/2 < r < 2 by proving the asymptotical optimality of order 2 digital nets for (1.1). On the 


other hand we would like to extend (1.3) to periodic Nikol’skij-Besov spaces with dominating 
mixed smoothness S^ q B(T d ), namely, 

QMC N (S^ q B(T d )) x IV r (logIV) (d_1)(1 ” 1/9) # N >2, (1.4) 

for 1/p < r < 2, see Definition |2.3| below. An immediate feature of these error bounds is the 
fact that the log-term disappears in case <7 = 1. Besov regularity is the correct framework 
when it comes to integrands of the form 


f(x) = max{0, g(x)} 


x E 


(1.5) 


so-called kink functions, which often occur in mathematical finance, e.g. the pricing of a Eu¬ 
ropean call option, whose pay-off function possesses a kink at the strike price, see e.g. [22l 
Chapter 1], In general, one can not expect Sobolev regularity higher than r = 3/2. However, 
when considering Besov regularity we can achieve smoothness r = 2. Indeed, the simple ex¬ 
ample f(t ) = max{0,f — 1/2} belongs to -£^^([0,1]) while its Sobolev regularity H s is below 
s = 3/2. In a sense, one sacrifices integrability for gaining regularity. Looking at the bound 
G3 above, we see that cubature methods based on order 2 digital nets benefit from higher 
Besov regularity while the integrability p does not enter the picture. 

Apart from that, spaces of this type have a long history in the former Soviet Union, see 
03 GEL 03j 49j and the references therein. The scale of spaces S r pq B(T d ) contains several 
important special cases of spaces with mixed smoothness like Holder-Zygmund spaces (p = 
q = oo), the above mentioned Sobolev spaces (p = q = 2) and the classical Nikol’skij spaces 
(q = oo). Note that Sobolev spaces S r p H(T d ) with integrability 1 < p < oo and r > 0 
are not contained in the Besov scale. They represent special cases of Triebel-Lizorkin spaces 
S/ tq F(T d ) if q = 2. However, classical embedding theorems allow to reduce the question for 
Q,MC N (SpH(T d )) to (1.4) in the case of “large” smoothness r > max{l/p, 1/2}, see Corollary 


1 These spaces are sometimes also referred to as Korobov spaces. 
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5.6 below. For a complete study of asymptotical error bounds (including the case of “small” 


smoothness) of numerical integration in spaces SpH(T d ) and Sp q F(T d ) we refer to the recent 
preprint [53]. See also Remark 5.7 below. 


The by now classical research topic of numerically integrating periodic functions goes back 
to the work of Korobov |28| . Hlawka m, and Bakhvalov [3] in the 1960s and was continued 
later by Temlyakov, see Hamms], Dubinin 03 H8] and Skriganov [54]. See also the survey 
articles Temlyakov m and Novak m ■ In particular, Temlyakov ggmi used the classical 
Korobov lattice rules in order to obtain for 1/2 < r < 1 


as well as 


N~ r (logN)^ < QMC N (S r 2>00 B(T d )) < N~ r (log N )(*-i)(r+i/2) 


N~ r (log iV)( d_1 )/ 2 < QMC N (H r mix (T d )) < N-^logN)^-^ 


( 1 . 6 ) 


(1.7) 


for N > 2. In contrast to the quadrature of univariate functions, where equidistant point grids 
lead to optimal formulas, the multivariate problem is much more involved. In fact, the choice 
of proper sets X.y c of integration nodes in the h-dimensional unit cube is the essence of 
“discrepancy theory” ini 12] and connected with deep problems in number theory, already for 
d = 2. 

Recently, Triebel [ED EZ] and, independently, Dung [19j brought up the framework of 
tensor Faber bases for functions of the above type. The main feature is the fact that the basis 
coefficients are linear combinations of function values. The corresponding series expansion is 
thus extremely useful for sampling and integration issues. Triebel was actually the first who 
investigated cubature formulas for spaces ST B(Q d ) of functions on the unit cube Q d := [0, l] d . 


By using more general cubature formulas of type (1.9) below (with non-equal weights) and 
nodes from a sparse grid Triebel obtained the two-sided estimate 

IV-moglV)^-^ 1 - 1 /*) < Intjv (S r p , q B(Q d )) < N-^logN)^ 1 ^ 1 - 1 ^ (1.8) 

if 1 < p, q < oo and 1/p < r < 1 + 1/p. Here, Intjv denotes the optimal worst-case integration 
error where one admits general (not only QMC) cubature formulas of type 


N 


Ajv(A/v, /) := ^ A if (a 


(1.9) 


2=1 


In contrast to Sp q B(T d ), the space Sp q B{Q d ) consists of non-periodic functions on Qd ■ = 
[0,1/Z The questions remain how to close the gaps in the power of the logarithms in (1.6), 
<0>, and (|1.8|) and what (if existing) are optimal QMC algorithms? 


This question has partly been answered by the first and second named authors for a subclass 
of Sp q B(Qd ) with 1/p < r < 1, namely those functions Sp q B(Qdy with vanishing boundary 


values on the “upper and right” boundary faces, by showing that the lower bound in (1.8) 


is sharp for quasi-Monte Carlo methods based on Chen-Skriganov points, see mmmm- 
Furthermore, together with M. Ullrich the last named author recently observed, that the 
classical Frolov method is optimal in all (reasonable) spaces S” q B(Qd) n and Sp q F(Qd) n of 
functions with homogeneous boundary condition, see [53]. Note, that Frolov’s method is an 


equal-weight cubature formula of type (1.9) with nodes from a lattice (not a lattice rule). In 
a strict sense, Frolov’s method is not a QMC method since the weights A i do not sum up to 1. 

In this paper we investigate special QMC methods for periodic Nikol’skij-Besov spaces 
on and answer the above question partly. The picture is clear in case d = 2, i.e., for 


spaces on the 2-torus T 2 . In fact, we know that the lower bound in (1.8) is even sharp for all 
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r > 1/p, see [48L 56j [20]. The optimal QMC rule in case 1/p < r < 2 is based on Hammersley 
points [56], which provide the optimal discrepancy in this setting [23] . This paper can be 
seen as continuation of [56] for the higher-dimensional situation by adopting methods from 

[Ml EH EH [33]- 

We will prove the optimality of QMC methods based on order 2 digital nets in the framework 
of periodic Besov spaces with dominating mixed smoothness if 1/p < r < 2. Due to the 
piecewise linear building blocks, we can not expect to get beyond 1/p < r < 2 with our 
proof method even when taking higher order digital nets. Therefore, this restriction seems to 
be technical and may be overcome by using smoother basis atoms like piecewise polynomial 
B-splines [19]. 

We illustrate our theoretical results with numerical computations in the Hilbert space case 
// r 7 ni x(T d ) in several dimensions d and for different smoothness parameters r. In the case of 
integer smoothness we exploit an exact representation formula for the worst-case integration 
error of an arbitrary cubature rule. A numerical evaluation of this formula indicates that the 
results in Theorem 5.3 below keep valid for r = 2. The comparison with other widely used 
cubature rules such as sparse grids and Halton points in all dimensions and Fibonacci lattices 
in dimension d = 2 shows that order 2 digital nets perform very well not only asymptotically 
but already for a relatively small number of sample points. Finally, we consider a simple test 
function which is a tensor product of univariate functions of the form (1.5). Expressing the 
regularity of such functions in Besov spaces of mixed smoothness allows the correct prediction 
of the asymptotical rate of the numerical integration error which is verified by our numerical 
experiments. However, the applicability of order 2 nets to real-world problems from option 
pricing etc., where the kinks are not necessarily axis aligned, requires further research. 

The paper is organized as follows. In Section 2 we introduce the function spaces of interest 
and provide the necessary characterizations and properties. The classical definition by mixed 
iterated differences will turn out to be of crucial importance. In Section 3 we deal with the 
Faber and Haar basis, especially with their hyperbolic (anisotropic) tensor product. The main 


tools represent Propositions |3.4| and 3.5 where the function space norm is related to the Faber 
coefficient sequence space norm and vice versa. In Section 4 we recall Dick’s construction 
of higher order digital nets and compute the Haar coefficients of the associated discrepancy 
function. We continue in Section 5 by interpreting the Haar coefficients of the discrepancy 
function in terms of integration errors for tensorized Faber hat functions. Combining those 
estimates with the Faber basis expansion and the characterization from Section 3 we obtain 
our main results in Theorem 5.3. Finally, Section 6 provides the numerical results. 

Notation. As usual N denotes the natural numbers, No = N U {0}, N_i = No U {—1}, Z 
denotes the integers, M the real numbers, and C the complex numbers. The letter d is always 
reserved for the underlying dimension in M rf , Z d etc. We denote by ( x, y) the usual Euclidean 
inner product and inner products in general. For a £ M we denote a+ := max{a, 0}. For 
0 < p < oo and x G we denote \x\ p = (Y/i=i \%i\ p ) 1 ^ v with the usual modification in the case 
p = oo. We further denote x + := ((xi) + ,..., (xd)+) and |x|+ := |x + |i. By {x\,... ,Xd) > 0 we 
mean that each coordinate is positive. By T we denote the torus represented by the interval 
[0,1], where the end points are identified. If X and Y are two (quasi-)normed spaces, the 
(quasi-)norm of an element x in X will be denoted by ||x|X||. The symbol X Y indicates 
that the identity operator is continuous. For two sequences a n and b n we will write a n < b n if 
there exists a constant c > 0 such that a n < cb n for all n. We will write a n x b n if a n < b n 
and b n < a n . 
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2 Periodic Besov spaces with dominating mixed smoothness 

Let T d denote the d-torus, represented in the Euclidean space by the cube T d = [0, l] d , where 
opposite points are identified. That means x, y £ are identified if and only if x — y = k, 
where k = (k\, ..., kd) £ Z d . The computation of the Fourier coefficients of an integrable 
d-variate periodic function is performed by the formula 

f(k) = f f{x)e~ i2nk - x dx , ke Z d . 

J T d 

Let further denote L p (T d ), 0 < p < oo, the space of all measurable functions / : T d —> C 
satisfying 

\\f\\p = {J Td i/( x )i Pdx ) /p < 00 

with the usual modification in case p = oo. The space C'(T d ) is often used as a replacement for 
Zoo(T d ). It denotes the collection of all continuous and bounded periodic functions equipped 
with the Loo-topology. 

2.1 Definition and basic properties 

In this section we give the definition of Besov spaces with dominating mixed smoothness on 
T d based on a dyadic decomposition on the Fourier side. We closely follow [43l Chapter 2], 
To begin with, we recall the concept of a dyadic decomposition of unity. The space C“(M rf ) 
consists of all infinitely many times differentiable compactly supported functions. 

Definition 2.1. Let 3>(M) be the collection of all systems ip = {99 n (x)}^L 0 Cl satisfying 

(i) supp (po C {x : \x\ < 2} , 

(ii) supp tp n c {x : 2 n ~ 1 < |x| < 2" +1 } , n = 1,2,, 

(in) For all f E No it holds sup 2 nl \D e ip n (x)\ < q < oo , 

x,n 

oo 

(iv) Fn{x) = 1 for all i£l. 

71=0 

Remark 2.2. The class <f>(M) is not empty. We consider the following standard example. Let 
<po(x) E 5(M) be a smooth function with p>o(x) = 1 on [—1,1] and <po(x) = 0 if \x\ > 2. For 
n > 0 we define 

ip n (x) = vo(2~ n x) - po(2~ n+1 x). 

It is easy to verify that the system = {^ n (.T)}^T 0 satisfies (i) - (iv). 

Now we fix a system <p = {<£>n} ne z ^ where we put <p n = 0 if n < 0. For j = (j i,... ,jd ) E 

7L d let the building blocks fj be given by 

f j (x)=Y^Fh(ki)---Tj d (k d )f(k)e i2 ^ x , x &T d ,j E . (2.1) 

kez d 
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Definition 2.3. (Mixed periodic Besov and Sobolev space) 

(i) Let 0 < p, q < oo and r > a p := ( 1/p — 1)+. Then Sl p q B(T d ) is defined as the collection of 
all f G Li(T d ) such that 

1/9 


\\f\S r m B(T d W := 


£ 2b1in7 ll/ 

jeNg 


* 11 ? 


( 2 . 2 ) 


is finite (usual modification in case q = oo). 

(ii) Let 1 < p < oo and r > 0. Then SpH(T d ) is defined as the collection of all f G L p (T rf ) 
such that 

\\f\Sr p H(T d )\r := ||( X) 2 |j ' l|r2 |/ J -(*)| 2 ) 1/2 

jSNg 

is finite. 

Recall, that this definition is independent of the chosen system ip in the sense of equivalent 
(quasi-)norms. Moreover, in case rninjjj, q} > 1 the defined spaces are Banach spaces, whereas 
they are quasi-Banach spaces in case min{p, q} < 1. For details confer | [43l 2.2.4]. In this paper 
we are mainly concerned with spaces providing sufficiently large smoothness (r > 1/p) such 
that the elements (equivalence classes) in S r p q B(T d ) contain a continuous representative. We 
have the following elementary embeddings, see [431 2.2.3]. 

Lemma 2.4. Let 0 < p < oo, r G M, and 0 < q < oo. 

(i) If e > 0 and 0 < v < oo then 

S r + q e B(T d ) ^ Sp jV B(T d ). 

(ii) If p < u < oo and r — 1/p = t — 1/u then 

S r p , q B(T d ) ^ Sl q B(T d ). 

(in) If r > 1/p then 

S r Ptq B(Y d ) C{T d ). 

(iv) If 1 < p < oo and r > 0 then 

^p,min{p,2}^C^b -+ S r r H (T J ) -4 S’^^B . 

(v) If 2 < p < oo and r > 0 then 

S r p H(T d ) S r 2 H(T d ) = H r mix (T d ) = SUB{ T d ). 


2.2 Characterization by mixed differences 


In this subsection we will provide the classical characterization by mixed iterated differences 
as it is used for instance in |fj. The main issue will be the equivalence of both approaches, 
the Fourier analytical approach in Definition |2.3| and the difference approach, see Lemma |2.7| 
below. We will need some tools from Harmonic Analysis, the Peetre maximal function and 
the associated maximal inequality, see [43^ 1.6.4, 3.3.5]. For a > 0 and b = (&i,..., bd) > 0 we 
define the Peetre maximal function Pb <a f for a trigonometric polynomial /, i.e., 


Pb,af{x) := sup 


(1 + bi\x\ 


\f(y)\ _ 

Vi\) a •■■(! + bd\x d - y d \) a 


The following maximal inequality for multivariate trigonometric polynomials with frequencies 
in the rectangle Qb := [—&i, 6i] x ... x [— b d , b d \ will be of crucial importance. 
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Lemma 2.5. Let 0 < p < oo, 6 = (b\, ..., 6^) > 0, and a > 1/p. Let further 

f 

\ki\<bi 

be a trigonometric polynomial with frequencies in the rectangle Q b . Then a constant c > 0 
independent of f and b exists such that 


E 


f(k)e 


2'Kik-x 


\\P b ,af\\p<c\\f\\ p . 

Now we introduce the basic concepts of iterated differences A f/(f,x) of a function /. For 
univariate functions / : T — > C the mth difference operator A™ is defined by 

m , v 

A^(/,x):=^(-ir^( m )f(x+jh) , x£T,fte[0,l]. 

The following Lemma states an important relation between mth order differences and 
Peetre maximal functions of trigonometric polynomials, see [5H, Lemma 3.3.1]. 

Lemma 2.6. Let a, b > 0 and 

f = J2 f(k)e 2 * ikx , 1 ST, 

| fc |<6 

be a univariate trigonometric polynomial with frequencies in [—6,6]. Then there exists a con¬ 
stant c > 0 independent of f such that for every h £ M 

|A™(/,x)| < cmin{l, |6/i| m } max{l, \bh\ a }P bya f{x) , x£T. (2.3) 


In order to characterize multivariate functions we need the concept of mixed differences 
with respect to coordinate directions. Let e be any subset of {l,...,d}. For multivariate 
functions / : T cl —> C and h £ [0, l] rf the mixed (m, e)th difference operator A)“’ e is defined by 

n A E.i “ d Ar*=id, 

zGe 


where Id / = / and A™ ■ is the univariate operator applied to the i-tli coordinate of / with 
the other variables kept fixed. Let us further define the mixed (m, e)th modulus of continuity 
by 


w m(/^)p : = SU P I 

\hi\<ti,i£e 


An/, on* 


t€ [0,l] c 


(2.4) 


for / £ L p (T d ) (in particular, w^(/, t) p = ||/|| p ). We aim at an equivalent characterization 
of the Besov spaces Sf q B(T d ). The following lemma answers this question partially. There 
are still some open questions around this topic, see for instance [43 . 2.3.4, Remark 2]. The 
following Lemma is a straight-forward modification of [54i , Theorem 4.6.1]. 


Lemma 2.7. Let 1 < p < oo, 0 < q < oo and m £ N with m > r > 0. Then 


where 


||/|^ <?J B(T d )rx||/|^ (?J B(T d )||( m ) , / £ Li(T d ), 

\\f\S r m B(T d )\\^ := [ Y, 2 rljliq u< j \f,2~X ll/q 

jeNg 


(2.5) 


In case q = oo the sum in (2.5) is replaced by the sup over j. Here, e(j) = {i : j, / 0}. 


7 




Proof. This assertion is a modified version of [5T, Theorem 4.6.2] for the bivariate setting. 
Let us recall some basic steps in the proof. The relation 


ll/l < s; gJ B(T rf )||( m )<c 1 ||/|5; i9J B(T d )r 


is obtained by applying m Lemma 3.3.2] to the building blocks fj in ( |2.1[ ), which are indeed 
trigonometric polynomials, and using the proof technique in |5l, Theorem 3.8.1]. 

To obtain the converse relation 

l\f\S r p , q B{T d )r <C 2 \\f\Sl .£(¥*) ||W 


we take into account the characterization via rectangle means given in Theorem 4.5.1]. 
We apply the techniques in Proposition 3.6.1 to switch from rectangle means to moduli of 
smoothness by following the arguments in the proof of Theorem 3.8.2. It remains to discretize 
the outer integral (with respect to the step length of the differences) in order to replace it by 
a sum. This is done by standard arguments. Thus, we almost arrived at (2.5). Indeed, the 


final step is to get rid of those summands where the summation index is negative. But this 
is trivially done by omitting the corresponding difference (translation invariance of L p ) such 
that the respective sum is just a converging geometric series (recall that r > 0). ■ 


Remark 2.8. By replacing the moduli of continuity (2.4) by more regular variants like integral 


means of differences we can extend the characterization in Lemma |i?.7| to all 0 < p < oo 
and r > (1/p — 1) + , see also Remark 3.6 below. 


3 Haar and Faber bases 



j e {o, 1} 


Figure 1: Univariate hierarchical Faber basis on T for levels j £ {0,1} and their union. 


3.1 The tensor Haar basis 

Let h(x) denote a piecewise constant step function on the real line given by 


h(x ) 


1 : 0 < x < 1/2, 

-1 : 1/2 < x < 1, 

0 : otherwise. 


For j 6 No and k 6 B j := {0,1,..., 2 J — 1} we put 

hj, k ■= • ~k). 
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Clearly hj >k is supported in [0,1] for j G No, k G B j . Let now j = (j \,..., jd) G Ng and 
k G Dj := B^ x ... x D j d . We denote by 

hj,k(x i,...,x d ) := h jlikl {x{)---h jdjkd {x d ) , (x 1 ,...,x d ) G [0, l] fZ , 

the tensor Haar function with respect to the level j G Ng and the translation k G and 

t J ‘j,k(f)= / f(x)h jk (x)dx 
J[o,l] d 

the corresponding Haar coefficient for / G Li([0, l] d ). 

3.2 The univariate Faber basis 

Recently, Triebel 15U (52] and, independently, Dung [T9] observed the potential of the Faber 
basis for the approximation and integration of functions with dominating mixed smoothness. 
The latter reference is even more general and uses so-called B-spline representations of func¬ 
tions, where the Faber system is a special case. We note that the Faber basis also plays 
an important role in the construction of sparse grids which go back to [4S] and are utilized in 
many applications for the discretization and approximation of function spaces with dominating 
mixed smoothness, see e.g. m esi 

Let us briefly recall the basic facts about the Faber basis taken from [511 3.2.1, 3.2.2], 
Definition 3.1. Let v(x) be the Loo-normalized, integrated Haar function h, i.e., 

v{x) := 2 f h(t) dt , iGM, (3-1) 

Jo 

and for j G No, k G Dj 

Vj,k = v(2 j ■ -k ). (3.2) 

For notational reasons we let U-po := 1 and obtain the Faber system 

F ■■= {v jtk : j G N_i, k G Bj} . 

Faber [21] observed that every continuous (non-periodic) function / on [0,1] can be repre¬ 
sented as 

1 oo 2-7 — 1 

f(x) = /(0) • (1 - x) + /(1) • x - g A 2-7 -i (/> 2_J k ) v j,k(x) (3-3) 

j= 0 k= 0 

with uniform convergence, see e.g. m Theorem 2.1, Step 4], Consequently, every periodic 
function from C'(T) can be decomposed via the (1-periodic extended) system F such that 

1 OO 2-7 -1 

/(*) = /(°) - 2 A 2-7-i (/. 2 ~ J k)v^k{x) (3.4) 

j=0 k= 0 


with convergence in C(T). 
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3.3 The tensor Faber basis 


Let now f(x i,... ,xa) be a d-variate function / G C(T d ). By fixing all variables except X{ we 
obtain by g{- ) = f(x\, ..., i, •, Xj+ 1 ,..., x^) a univariate periodic continuous function. By 
applying (3.4) in every such component we obtain the representation in C{T d ) 


/(*)= d lk(f) y j,k( x ) > xGT d , 

jeNl 1 fceBj 


(3.5) 


where 

^(xi, ■ • • ,x d ) := v jljkl (xi) ■ ■ -v jdtkd {x d ) , j G G Dj , 

and 

4 fe (/) = (-2)- |e W | A^g! 1) (/,^fc) , jeNi^eD,, (3.6) 

Here we put e(j) = {i : j* 7^—1} and Xj^ = ( 2 _ ^ 1 )+^ 1; ... ; 2 ~^ d ' )+ kd) ■ 

Our next goal is to discretize the spaces S r pq B(T d ) using the Faber system {vj^ : j G 
Nl l5 k G Dj}. We obtain a sequence space isomorphism performed by the coefficient mapping 
d 2 j k (f) above. In [51, 3.2.3, 3.2.4] and |19l Theorem 4.1] this was done for the non-periodic 
setting S pq B(Qfi). Our proof is completely different and uses only classical tools. This makes 
the proof a bit more transparent and self-contained. With these tools we show that one 
direction of the equivalence relation can be extended to 1/p < r < 2. 


Definition 3.2. Let 0 < p, q < oo and r G M. Then s pq b is the collection of all sequences 
Nl-pfceBj su °h that 


y,k I 


'■ii - [ E 2b‘h P-i/FL 

je N£j 


ke ID),- 


q/p- 


1/9 


is finite. 

Lemma 3.3. Let 0 < p, q < oo and r G M. The space s pq b is a Banach space if min{p. q } > 1. 
In case min{p, q} < 1 the space s p q b is a quasi-Banach space. Moreover, if u := min{p, q, 1} 
it is a u-Banach space, i.e., 


iiA+pi4^r<iiAi4 g 6r + iH4 q c 


A, h £ s p,qb ■ 


Proposition 3.4. Let 1/2 < p < oo, 0 < q < oo and 1/p < r < 2. Then there exists a 
constant c > 0 such that 

IKy/K^II < c||/|s;,,B(r')r (3.7) 

for all f € C(T d ). 


Proof. The main idea is the same as in the proof of Lemma 2.7 We make use of the decompo¬ 
sition (2.1) in a slightly modified way. Let us first assume 1 <p,q< oo. The modifications in 
case minjp, q} < 1 are straight-forward. For fixed j G we write / = fj+t- Putting 

this into (3.7) and using the triangle inequality yields 


E E d Uf)vi*w 

je. n£j 




91 1/9 
V- 


< 


E [ E 2rW,, || E 


i&z d je Nlj 


fcSB,- 


1/9 


(3.8) 
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Let us continue in deriving a point-wise upper bound for the absolute value of the function 


F jA x ) '■= Y d lk(fj+e) v j,k( x ) ■ 

keOj 


Clearly, we have 

\ F jA x )\ < \ d lk(fj+e)\ fz \ A 2 2 -uh)(fj+e, x j,k)\ ( 3 - 9 ) 

whenever x G 2~^ d> +{k\ + 1)] x ... x [2~^ d ^+kd,2~^ d ^+ (kd + 1)] • Let us estimate 

the iterated differences A l-u+i)(fj+i^ x j,k) one by one. For i G e(j) we have for x such that 
lx.,-, h — x\ < 2~ 3i the bound 

I Jii^i I — 


1^2 ~(h+ 1 '> Ai x ji,ki) | 


< sup \g(x + y)\ 
\v\<2~ j i 


< 


sup 

\y\<2~ j * 


\g{ x + y)\ 

(1 + 2L|y|)« 


< p 2H. a a( x ) 


for a univariate continuous function g. In case i ^ e(j) we have j t = — 1 and 


\g(0)\ < sup \g{x + y)\ < sup 
|y|<i M<i 


\a{x + y)\ 
(1 + 2%|)“ 


< p 2n, a a( x ) ■ 


if £ j > 0 then, by definition, 


p 2n, a a( x ) < 2 ' ap 2n+k^aA) ■ 


(3.10) 


On the other hand, the estimate in Lemma 2.6 gives in case i G e(j) for a univariate trigono¬ 
metric polynomial g ji+ At) = EfceZ <Pji+ti(k)g(k)e 2mkt 


\\-(ji+Aaji+in x ji,ki)\ ~ min{l,2 'jmaxll^ ia }P 2 ^ +ii a5ii+£ .(x i . ifci ) 
If ii < 0 and |x — < 2 _:,i this reduces to 

l 2 ^2 _ 0’i+ 1 ) (3ji+£ii x ji,ki)\ 2 1 P 2 e i+ii,a9ji+£i ( x ) • 


(3.11) 


Note, that in case i ^ e(j) there is nothing to prove since gj i +e i = 0. Applying the point-wise 
estimates (3.10) and (3.11) to the right-hand side of (3.9) we obtain 

\ p jA x )\ < p 2e+j,afj+e(x) n min{2 2/; ' : , 1} max{2 £i “, 1} , 

iGe(j) 


where 2 i+3 := (2^ 1+Jl ) ..., 2 ^ d+3d ). Using the Peetre maximal inequality, Lemma 2.5 yields 

(3.12) 


\ p jA\p ^ \\ p 2t+i,afj+e\\ P ■ II min (2 2 'M} max{2^ a , 1} 

■i£e(j) 

\\fj+i\\p II min{2 % ,l}max{2^ a ,l}, 
iSe(i) 


< 

r^j 


whenever a > 1/p. If r > 1/p we can choose 

l/p<a<r<2. 

Therefore, if i G 


£ 2 rUI,, iiF«ii? s E ^"li/^elK. 


(3.13) 


(3.14) 


2=1 
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where for ngZ 


A n = 


2 ( 2 —r)n 
2 (a—r)n 


n < 0, 

n > 0. 


(3.15) 


Under the condition ( 3. 13| ) it follows from (3.15) that there is a 5 > 0 such that A n < 2 <5 l n l 
and hence 

E 2rljhq \\ F jA g p < 2-^' 1 ||/| S r m B(T d )\\i. 


Plugging this into ( |3.8[ ) yields ( |3.7[ ) . ■ 

Let us prove the converse statement. The version below slightly differs from its 2-dimensional 
counterpart given in [56| although the proof technique is the same. We observed that the re¬ 
striction r > 1/p is actually not required. 

Proposition 3.5. Let 1 < p < 00 , 0 < q < oo, 0 < r < 1 + 1/p. Then there exists a constant 
c > 0 such that 

ii/is;,,B(T < )ir < 0.16) 

for all f € C(T d ) with finite right-hand side (3.16) . 


Proof. We use the characterization in Lemma |2.7[ which says that 

< ll/|S; s B(T , )ll <m) 

for some fixed m > 2. Let us assume 1 < q < 00 . The modifications in case q < 1 are 
straight-forward. Similar as done in the previous proof we obtain by triangle inequality 


jeNg 


1 1/9 


< 




1 i/q 


(3.17) 


je Ng 


where we put (in contrast to the previous proof) 


fj( x ) = E d lk(f) v jA x ) i 

fceB,' 


with fj = 0 if j (f . We exploit the piecewise linearity of the basis functions Vj + e^ together 
with the at least second order differences in In fact, let us consider the 

variable xi. For i\ < 0 and \h\\ < 2~ n the difference (vj 1 +^ 1) fc 1 ,xi) vanishes unless xi 
belongs to one of the intervals Il,Im,Ir given by II '■= {x E T : |x — 2 ~^ 1 +h)k l \ < |hi|}, 
I M :={x& T : |x - 2 ~(n+h)( kl + 1/2 )| < |/n|}, an d Ir ■- { x G T : |x - 2~^ + A){k l + 1)| < 
|hi|}. Furthermore, in case t\ < 0 it is easy to verify that 

l A Ai( u ii+*i,fci.*l)l < 2* , xi € II U Im U Ir . 

In particular, as a consequence of \I P U Im U Ir\ < \h\\ < 2 _jfl we obtain 

[ |(A| 1 ^ 1+ ,„ ll )(i 1 )|>’dx 1 <2^2-* (3.18) 

J T 

in case l\ < 0. In case i\ > 0 we use the trivial estimate 

/ KA^+^KmP-dx! < . (3.19) 

J T 
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Now we combine the component-wise estimates in (3.18) and (3.19) to estimate \-e, 2 ] ) p 

from above. Indeed, using the perfect localization property of the basis functions we obtain 


< (2- fa+,) y i^ + „(/)r) 1/, 'n{ 

k£Bj +e i =1 ^ 


2 b( 1 + 1 /p) . £ i< Q i 


li> 0 • 


Now, similar as in (3.14) in the previous proof we see for t E IT 


jeNg 

< ^ 2 (r " 1 / p )l i+ ^ 19 


ieNg 

which results in 




?/p f 2 £i ( 1+1/,p_? ')' ? : t 

2 ~riiq 


< 0 , 

i>0, 


(3.20) 


E 2 rljlirq u< j Hf j+ e,2-x < 2-^i^ii4 fe (/)i^6r, 

jeNg 


where we used that 0<r<l + l/p. Plugging this into (3.17) concludes the proof. 


Remark 3.6. The restriction p > 1 in Proposition 3.5 can be removed. Note, that this 
restriction is caused by the difference characterization in Lemma [¥7?\ which can be extended to 
0 < p, q < oo and m > r > (1/p — 1) + , see [5/1 Theorem 4-5.1], by using rectangle means of 
differences, 


K e m (f,t) p := 


IA 






(3.21) 


instead of the mixed moduli of continuity (2.4). In other words, if 

0 < p, q < oo and ( 1/p — 1) + < r < 1 + 1/p 
it holds with m . > max{2,1 + 1/p} 

\\f\S r m B(T d )r < \\f\S r Ptq B(T d )\\P < \\dl k (f)\s r m b\\ 


(3.22) 


for all f E C(T d ) with finite discrete quasi-norm Wd'j k (f)\s p q b\\ . The quasi-norm in the middle 
of (3.22) represents the counterpart of (2.5), where (2.4) is replaced by (3.21) . Note, that the 
restriction r < 2 does not occur here in case p <1. 


4 Discrepancy of order 2 digital nets 
4.1 Digital (t,n,d)~ nets 

We quote from (TJ Section 4] to describe the construction method of order a digital (t, n, d)-nets 
which in case a = 1 are original digital nets from [36] but in this form they were introduced in 

0. 

For s, n E N with s > n let C ±,..., Cd be s x n matrices over F 2 . For v E {0,1,..., 2 n — 1} 
with the dyadic expansion n = vq + n±2 + ... + z/ n _i2 n_1 with digits vq, Vi, ..., v n -i E {0,1} 
the dyadic digit vector u is given as v = (i/q, v \,..., v n - i) T £ F£. Then we compute C{i> = 
( x i,v.i , x i,y, 2 , ■ ■ •, a^,s) T G F| for 1 < i < d. Finally we define 

%i,v = T .V.2 2 T • • • T 2 E [0,1) 


13 
















and x v = (x \, v ,..., Xd,v)- We call the point set V n = {£o> xi, • • ■, £ 2 "-i} a digital net over F 2 . 

Now let a G N and suppose s > crn. Let 0 < t < an be an integer. For every 1 < i < d 
we write Ci = (c^i,..., Ci >s ) T where c^i,..., Cj )S G Fl> are the row vectors of Q. If for all 
1 < Ai.i < • • • < A i tVi < s, 1 < i < d with 

Ai,i + • • • + + ... + A^i + ... + A c ; )m i n { 77d)(T } < an — t 

the vectors ci,a m , • • •, • 1 c LA d ,i, • ■ • > c d,\ d , Vd are linearly independent over F 2 , then V n 

is called an order a digital (t,n,d )-net over F 2 . 

The quality parameter t and the order a qualify the structure of the point set, the lower t 
and the higher a - the more structure do the point sets have. 

Lemma 4.1. Let V n be an order 1 digital ( t , n, 6)-net then every dyadic interval of order n — t 
contains exactly 2 4 points ofV n . 

Therefore (t, n, ri)-nets are also (t + 1, n, ri)-nets and order a + 1 nets are also order a nets 
(with even lower quality parameter). In particular every point set V n constructed with the 
method described above is at least an order a digital (an, n, d)-net. We refer to [6j and [7] for 
more details on such relations. 

We need the following fact concerning projections of digital nets (0 Theorem 2]). 

Lemma 4.2. Let V n be an order 2 digital (t,n, d)-net. Let further Ii := {ii, ■ ■ ■ ,ie} C 
{1,... ,d} be a fixed set of coordinates. Then the projection V n (Ii) C [0, Vf of the set V n 
on the coordinates in I? is an order 2 digital (t,n,i)-net. 

Now we quote explicit constructions of higher order digital nets. We will only briefly 
describe the method, for details consult m,m and [8j. The starting point are order 1 digital 
(t',n,ad )-nets and the so called digit interlacing composition 

■■ [ 0 , ir ->• [ 0 , 1 ] 

(XI, • • ■ , X a ) ^ f; tr,a2- r - {a ~ 1)a , (41) 

a =1 r =1 

where £r,i, f r , 2 - ■ ■ ■ are the digits of the dyadic decomposition of x r . The digit interlacing is 
applied component wise on vectors, namely 

(x\, , Xfyd) 1 ^ (^u (xi , ••• , Xaj, , & ( ?(X( c i—l)(T+l j •.. j yd) ) • 

element wise. Suppose that V n is an order 1 digital (t',n,ad)- net. Then S>^(T) is an order a 

digital (t,n,d )-net with t = at' + da(a — l)/2. Therefore, it is possible to construct order 2 
digital ( t , n, el)-nets. 

4.2 The discrepancy function and its Haar coefficients 

Let N be a positive integer and let V be a point set in [0, l] d with N points. Then the 
discrepancy function D-p is defined as 

D v( x ) = jj^2x{o,x]( z ) - x i--' x d (4.2) 

z&V 

for any x = (xi ,..., Xd) G [0, l] d . By X[o,x] we denote the characteristic function of the 

interval [0,x] = [0,aq] x ... x [0,Xrf], so the term ^2 z X[o,x]( z ) is equal to jf(V n [0,a;]). Dp 
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measures the deviation of the number of points of V in [0, x\ from the fair number of points 
7V|[0,x]| = N xi • • -x d . 

For further studies of the discrepancy function we refer to the monographs mmmm 
and surveys mmm- 

We will use Haar coefficients of the discrepancy function which are given by [331 Proposition 
5.7], 

Proposition 4.3. Let V n be an order 2 digital ( t,n,d)-net over ¥ 2 - Let further j 6 Nq and 
m 6 Bj. Then there exists a constant c = c(d) > 0 that satisfies the following properties. 

(i) If |j|i >n - \t/2] then 

H,m(DvJ < c2-^~ n+t / 2 

and /J‘j,m(D , p n ) < c2~ 2 ^' 1 for all but 2 n values of m . 

(ii) If |j'|i < n - ft/2] then 

Hm(D Vn ) < c2~ 2n+t (2n - t - 2\j\ l ) d ~ l . 

Remark 4.4. (i) It is in this proposition that the higher order property of the digital nets is 
needed. For a usual order 1 digital ( t,n,d)-net the main factor in the second estimate would 
only be 2~^^ 1 ~ n ~ t instead of 2~ 2n+t which is not sufficient to yield the right order of conver¬ 
gence in our results. 

(ii) In f24 y the author computed Haar coefficients of the dicrepancy function with respect to 
the two-dimensional Hammersley point set. In contrast to the current paper the author con¬ 
sidered the non-periodic situation Sf q B([(), l] rf ) and therefore had to deal with the cases where 
components ji of j equal —1 as well. This corresponds to the inner product of Dp n with the 
characteristic function X[o,i] the respective component i, see [2fl Theorem 3.1, (iv)]. This 
upper bound is essentially sharp, as shown in [2f\ Lemma 3.7], and responsible for the fact, 
that the results in \2Ifl can not be extended to r < 0. A d-dimensional counterpart for Chen- 
Skriganov points can be found in W- 


5 QMC integration for periodic mixed Besov spaces 


In the sequel we consider quasi-Monte Carlo integration methods for approximating the integral 
1(f) := fj d f(x) dx of a d -variate continuous function / G C'(T d ). More precisely, for a discrete 
set V C [0, l] d of N points we compute 


lN(V,f):=^J2f( Z ) > /^> 

zer 


where Fj denotes a class of functions from C(T a! ). 
ate Faber expansion (3.5) converges in C(T d ). We 
In( fP 1 f) ~ 1(f) ■ In fact, 


Assume that for / G Fj the multivari- 
consider the integration error i?jv(/) := 


\Rn(I)\ 


zGV 



f(x)dx 


E E d lm(f) jyr E V j,m( z ) 

j'gNl x mSBj zGV 

E J2 d UU)Cj,m(r) |, 


E E f d Vj,m(x)dx 

j€ Nl x rnGOj ^ 


(5.1) 
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where 


(5.2) 


Cj,m{V) :=- V Vj,m(z) - / Vj t m(x) dx , j £ N^, 771 £ Dj . 

iV * ' /Td 

z&P 1/1 

Let us first take a look at the second summand. 

Lemma 5.1. Let j £ and m £ Dj then 

[ Vj, m (x) dx = 2 _|j+l11 . 

JT d 

Proof. We use the tensor product structure of the Vj >m to compute 


h d 


u j,m 


d n d 

(x) d X = T[ v jumi ( Xi ) d Xi = TT 2 

i=l JT i =1 


b'i+i| — 2 _ lj , '+ 1 l 1 


The next Lemma connects the Haar coefficients fj,j >m (D-p) of the discrepancy function Dp 
with the numbers Cj tm (V) . 

Lemma 5.2. Let V C [0, l] d with ffV = N. 

(i) If j £ Nq and m £ Dj we have 

Hm(Dp) = (-1 ) d 2-\ih Cj , m (l>) , 

(ii) If j £ Nl 1 \ Ng and m £ Dj we have 

»j,m(Dp) = (-l) s 2-W+ Cj , m (V ), 


where V denotes the projection of V onto those s coordinates zi where ji —1. Moreover 
hj, m {Dp) is the Haar coefficient with respect to the s-variate function Dp, 

Proof. Let j £ Nq and m £ Bj. We compute pLj^ m {Dp). This involves two parts. We first 
deal with 


'[o,i] 


1 

d N 


Y X[0,x)(z)hj,m(x) dx = ^ Y [ X[z,l)( x ) h j,m( x ) dx 
z&V 2e p-'[0,l] d 


W] d 

d r l 


h Y II [ bjumtiy) dy. 


(5.3) 


z€V i= l J z i 


Let us deal with the univariate integrals on the right-hand side of (5.3). Clearly, for any 
i = 1, 


.,... ,d, 


-1 rzi rzi 

hji,mi (y) d y = - h H . m% (y) dy = - h(2 h y - mf) dy 

Zi JO Jo 

r2 J i Zi—rrii 

= —2~ Ji / h(r) dr 

Jo 


(5.4) 


This together with (|5.3|) yields 

[ I 

1 id N 




=-2 ji v(2 ji Zi - mi) = -2 h v jumi (zi) 


5 ^X[ 0 ,*)W^i,m(*) d ® = (-!) d 2 Y v l™( z ) ’ ze[0,l} d . 


z£P 


(5.5) 


z£V 
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It remains to compute 


d r i 


/ x 1 ■ ■ ■ x d hj, m (x)dx = TT / y h jumi {y) dy . 

■Ao,i] d f=1 Jo 


(5.6) 


Integration by parts together with (5.4) yields for i = 1,..., d 


"1 rl 

y hj^rm (y) d y = -2~ ]i / v jumi (y) d y 
) Jo 


which, together with (5.6), implies 


[ xi-■ ■ x d hj, m (x)dx = (-l) d 2 [ v j:m (x)dx (5.7) 

J[o,i]d 7[o,i] d 


Combining, (5.3), (5.4), and (5.7) yields the result in (i). The result in (ii) is a simple conse¬ 
quence of (i). ■ 

The following result represents our main theorem. 

Theorem 5.3. Let V n be an order 2 digital ( t , n, d)-net over F 2 . Then for 1 < p, q < 00 and 
1/p < r < 2 there exists a constant c = c(p, q , r,d) > 0 and we have with N = 2 n 

QMC N (S r pq B(T d ))< sup \I(f)-I N (V n ,f)\<c2 rt / 2 N- r (\ogN^ d - 1 ^ 1 - 1 ^ , nG N. 


Proof. Let / 6 Sp q B(T d ). By the embedding result in Lemma 2.4 /(ii),(iii) we see that 
/ G 1 B(T d ) for an e > 0. As a consequence of Proposition 3.4 we obtain that (3.5) 
converges to / in C(T d ) and therefore in L p (T d ). Then, by (5.1) together with twice Holder’s 
inequality we obtain 

\RnU)\ < E E l4m(/)Ci,m(Pn)| 


< 


[ E E K>.(/)i" 


d 


[E 


mgDj 




q/p- 


1/9 


^ ) \Cj } m(fPn 


t\q'/P 




w 


(5.8) 


me I 


ii/i^,s(t‘)ii • [ e 2- ir - 1/rm !•'( e 




m£D 7 - 


where we used Proposition |3.4| in the last step. In order to prove the error bound it remains 

q' /p' 


to estimate the second factor. Let us deal with 


bii<n-p/ 2 ] mdOj 

IS N d 


I 
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first. By Lemma 5.2 together with Proposition 4.3,(ii) we obtain 


x; 2-(r- i /p)bi i *'( ^ \ Cj ^y ,/p ' 

\j\l<n-\t/2\ meBj 

j 6Ng 

< ^ 2 — I- 7 1 1 (’— 2 ( —2ri +*)2 lb 1 1 2 1 1 / p/ ( 2 n — t — 2|j | i)( rf ~ 1 )'5 ,/ 

\j\l<n- ft/2] 

jeNg 

< 2 (- 2 n+t)g' ^2 2 - b1iP- 2 )9'(2n - t - 2|j|i) (d “ 1)9 ' 

|/|i<n-ft/2] 

ieNg 

n— ft/2] —1 

x2 (-2n+ty ^- 1 (2n-t-2£) (d_1) ' ?, 2" £(r " 2) ' ? '. 

e=o 


Putting M \= n — [~f/2] we obtain 


(5.9) 


,, , M—l 

2 -(r-l/p)|/| l9 '^ ^ | c . m |Py P x 2 (-2n+t)g' ^ t d ~ X (Af - ^(<*-l)9 , 2 -<(r-2)g / 
|/|l<n— ft/2] mgDj £=0 

jeNg 


^ 2 (—2n+t)g' 2 (2— r)Mq '1 


M—l 

x 2^- 2 ~ r ^ i ~ M ^ q ' (£/M) d ~ 1 (M — . 

t=o 


At this point we need the assumption r < 2 in order to estimate 


M—l 

Y 2( 2_r )(^ _M )9'(^/M) d_1 

^=0 


(Af - f) (,i ~ 1)9 ' < ^ 2 ^ (2 " r)9 '£ (d " 1)9 ' < C < 
n=o 


oo. 


This gives 


2 -(r-i/p)\j\iq'\c jm \P'Y /P < 2 {- 2n + t ) q '2^ 2 ~ r "> Mq 'M d ^ x x 2~ rnq ' 2 rq,t / 2 

b'|i<jr— ft/2] 

jeNg 


mGD 7 - 


Let us now deal with X)|j|i>n-ft/ 2 ] • By Proposition 4.1, (ii) and Lemma 5.2 we get 

Y 2 -b- 1 /P)bli'^ Y \CjmV'^ q ' /P ' 

b]i>n-ft/2] mgB, 


o 


< Y^ 2~^ r ~ l / p ^ iq ' 

\j\l>n-\t/2\ 
je Ng 




,A<?7p' 


m£Aj 


+ Y^ 2" (r “ 1/p )bii ,? ' 

|j|i>n-r*/2l 
'o 


E 


u .?) m I 


IP 


A «7p' 


jeNg 




where Aj d enotes the set of indices rri G Dj where I J)m nV n / 0- Clearly #Aj < 2 n . By 
Lemma 5.1 we directly obtain \cj jTn \ < 2 _ b’+ 1 li if m g D - \ Aj, whereas by Lemma 5.2 and 
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Proposition |4.3[ (i), |cj, m | < 2 n+t / 2 if m E Aj. This gives 


yy 2 — — i /p)IjI i 9 / 

Ii|i>n-r*/2l 

beNg 


E 

m£Aj 



< 


yy 2-( r - 1 /p)bii9 , 2 n, ?7p , 2(- n + t / 2 )7 

|j|l>n- rt/2] 


< 2^~ n+t / 2 '> q '2 nq '^ P ' 


E 


2-(r- 1 /p)|i|iq' 


U|i>n-rt/ 2 l 

jSNg 


< 2( _n+t / 2 ^ , 2 Tl ' 5,, / p, 2 _ ^ T ’ _1//p ' )Tl ' J, n d_1 x 2~ rnq ' 2 rq,t / 2 n d_1 


(5.10) 


where we used r > 1/p in the last step. Furthermore, 

,„a7/p' 


yE 2 _ 7- 1 /p)lili<?' 

bii>n-r</2i 

beNg 


E I c j,rn I 


< 


yp 2 _ 7- 1 /p)lili'? , 2~hli9'2lhk7p' 

iiii>»-rt/2i 

beNg 


< y^ 2 _ I j, I 1 7~ 1 /p+ 1 - 1 /p')'?' < 2 ~ nrq ' n d ~ l 
b'|l>n- rt/21 


beNg 


Putting everything together yields 


[E : 

beNg 


) -(r-l/p)h'|ig' 






mGDo 


i/7 


< 2 -rn 2 r 7 2 n ( d - 1 )( 1 ” 1 /'?). 


(5.11) 


It remains to consider the sum XOjgjp ,\Ng • I 11 f ac L we decompose 

E - E E • 

feNi x eC{l,...,d} jeNl 1 (e) 


where Nl j (e 
with Lemma 


= {j E Nlj : ji E No if i E e, and ji = — 1 if i ^ e}. By Lemma 5.2, (ii) together 


4.2 


we can estimate Yljew 1 ,( e ) by means of (5.11) and obtain for any fixed e / 

,\7/p'i W 


yP 2-P- 1 /p)b'|i7 
beNg(e) 


E l C 7-(^n)| P ' 


< 2 _rn n^ e l _1 )( 1_1 /9) . 


(5.12) 


mGI 


Note, that in case e = 0 we obtain c/ i -lyoC^n) = 0. Finally, (5.11) and (5.12) together 
yield 

9 /P 1 !/9 < 2 _ rn n (d-l)(l-l/g) ^ 


[ y^ 2 - 7-i/p)b'li7 
which concludes the proof. 


E iwnr 

mGBo 


Remark 5.4. (i) IFe comment on the question of extending the result in Theorem 5.3 to non¬ 
periodic Besov-spaces of mixed smoothness Sf q B([0,l] d ). Similar as in [24l [30!/ we are able 
to prove a corresponding bound for non-periodic spaces if r < 1 using order 2 digital nets. By 
Remark f.f./ii) and the correspondence in Lemma 5.2 this result may not extend to r > 1 in 
general. Note, that the Hammersley points represent an order 2 digital net in d = 2. Therefore, 
the correct order of QMC N (Sf q B([ 0, l] rf )) in the non-periodic situation r > 1 is open. 
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(ii) In order to integrate non-periodic functions one would rather use a standard modification 
of the QMC-rule (which then is not a longer a QMC-rule) given by 


n n ^(4 


Q%(f) ■■= iN(x N ,i$f ) = ^ 

3 = 1 


= 1 


N 


... ,fi(x J d )), 


(5.13) 


where 

rplp . 


: f i—^ (n ^'0*4)) • fW x i), • • •, ip(xd)) , f e ii([0, l] d ), x = (xi,..., x d ) E [0,1] (/ , 


2=1 


and f)' is a univariate sufficiently smooth bump with 0) = ifi/T) = 0. IfT d : SJ S([0, l] d ) —» 
Sp q B(T d ) is a bounded mapping, see (18/ and the upcoming paper [35], then the cubature 
formula Q ^ has the same order of convergence on S p>q B([ 0, l] d ) as LjvpGv, •) on the periodic 
space Sp q B(T d ). 


The result in Theorem 5.3 is optimal. In fact, the following lower bound for general cubature 
rules has been shown in HZ! and independently with a different method in [20] . 

Theorem 5.5. Let 1 < p, q < oo and r > 1/p. Then we have 

Int N {S r pq B{Y 1 )) > jV-^loglV)^- 1 ^ 1 - 1 ^) , NYN. 

By embedding, see Lemma |2.4[ (iv), (v) we directly obtain the following bound for the 
classes S r p H{T d ). 

Corollary 5.6. Let V n be an order 2 digital (t,n,d)-net over fV Then for 1 < p < oo and 
max{l/p, 1/2} < r < 2 there exists a constant c = c(p, q, r, d) > 0 and we have with N = 2 n 


QMC N (S r H(T d )) < 


p-v- ,, - SU P 

||/|5-H(Trf)||<l 


\I(f)-I N (V n ,f)\ <c2 rt / 2 N- r {\ogN)^ d - 1 ^ 2 , n G N . 


Proof. If p > 2 we use the embedding Lemma 2.4, (v) together with Theorem 5.3 Note, that 


we need r > 1/2 here. If 1 < p < 2 we use Lemma 2.4, (v) together with Theorem 5.3 


Remark 5.7. The case 2 < p < oo and 1/p < r < 1/2 is not covered by Corollary 5.6. 
This situation is often referred to as the “Problem of small smoothness”. It is not known how 
digital nets (order 1 should be enough) behave in this situation. Temlyakov [Jfyl was the first 
who observed an interesting behavior of the asymptotical error for the Fibonacci cubature rule 
in the bivariate situation in spaces S r p H( T 2 ). Recently, in m this behavior has been also 
established for the Frolov method in the d-variate situation. In fact, for spaces SfH(T d ) n with 
support in the unit cube Qd it holds for 1 < p < oo and r > 1/p 


Int N {S r v H{Q d ) u ) < N~ 


(logJV) (d_1)(1_r) 

(log lV)( d-1 )/ 2 \/log log N 
(log N) (rf- 1 )/ 2 


p > 2 A 1/p < r < 1/2, 
p > 2 A r = 1/2, 
r > max{l/p, 1/2} . 


We strongly conjecture the same behavior for S p H(T d ) where classical digital (order 1) nets 
give the optimal upper bound. 
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Figure 2: Worst-case errors of order 2 digital nets in H1 nix (T d ) for smoothness r = 1. 


6 Numerical experiments 

In this section we use the theory of reproducing kernel Hilbert spaces (RKHS) to explicitly 
compute the worst-case error for particular constructions of order 2 digital nets based on 
Niederreiter-Xing sequences in the case of integer smoothness r £ {1,2}. Moreover, we give 
numerical examples for the case of fractional smoothness, i.e. r = 


6.1 Worst-case errors in H^ lix (T d ) 

Let us recall that the Besov space coincides with the classical tensor product Sobolev 

space ix (T d ) := H r (T) (g)... (g) H r ( T) of functions with mixed derivatives of order r bounded 
in L 2 (T d ). Since for r > 1/2 the space ffj,^(T^) is a Hilbert space which is embedded in 
C(T d ), it is well-known [2] that for a given choice of an inner product (v)flj. there exists a 
symmetric and positive definite kernel K : T d x T d — > M that reproduces point-evaluation, i.e., 
it holds f(x ) = (/(■), K(-, x))h^ x for all x E [0, l) d and / £ H^- x (T d ). Then one can use the 
well-known worst-case error formula to compute the quantities 


II Rn\(H, 


r \*i|2 
mix / 11 


= sup 
l!/l^ixll<i 


N 


I T d J T d 


f(x) dx — N' \if(x l ) 

i Z ' 

i =1 

N .... 

K (x, y) dxdy — 2^2 Xi / K(x\y)dy + XX Xi\jK(x\ 

i= 1 i =1 i=l 


( 6 . 1 ) 


N N 


X J 


( 6 . 2 ) 


explicitly, if a point set Vn = {x 1 ,..., x N } of integration nodes and associated integration 
weights Ai,..., Atv £ K are given. In order to have a simple closed-form representation of the 
kernel K we choose the inner product of the univariate Sobolev space H r ( T) to be 


(f,g)Hr(T) = /(O)ff(O) + l 2vrfc l 2r f(k)g{k). (6.3) 

fcez\{ o} 
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Figure 3: Worst-case errors of order 2 digital nets in ifF x (T rf ) for smoothness r = 1. 


The induced norm is given by 

OO 

ii/Hh' = i/(o)i 2 + X i“i 2r i/mi 2 

k£ Z\{0} 



f(x)dx 



(6.4) 


where the last equality only holds for r £ N. Then the reproducing kernel K\ : T x T —>• M of 
H r ( T) is given by [57] 


Ki t r(x, y) :=1 + ^2 \2nk\ 2r exp(27rife(x - y)) 

fee z\{o} 

OO 

=1 + 2 ^22 | 27t/c| — 2r cos(27t k(x — y)). 
k=i 


(6.5) 


If r G N the kernel can be written as 

(- 1++ 1 

Ki, r (x,y) = 1 + V , . B 2r (\x-y\), (6.6) 

(2r)! 

where E> 2 r : [0,1] —> M denotes the Bernoulli polynomial of degree 2r. Since is the 

tensor product of univariate Sobolev spaces, the reproducing kernel of H ^ nW (T d ) is given by 
the product of the univariate kernels, i.e. 


d 

Kd,r( x > y) = II R \r( x j,yj), X,y eY l 

3 = 1 

reproduces point evaluation in i7+ x (T rf ). 

As an example we employ order 2 digital nets that are based on Xing-Niederreiter se¬ 
quences m , which are known to yield smaller t -values than e.g. Sobol- or classical Niederreiter- 
sequences [12] . For the special case of rational places this construction was implemented by 
Pirsic 023, see also [13] . It is known [38] that one obtains a digital (t,n,d)- net from a digital 
(t, n,d — l)-sequence {x°,..., x 2 " -1 } by adding an equidistant coordinate, i.e. 

{(L4)jn,.... L4 -i)J«» i/2") : i = 0,..., 2" - 1}, (6.7) 
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Figure 4: Worst-case errors for smoothness r = 2 and various cubature methods. 




Figure 5: Worst-case errors for smoothness r = 2 and various cubature methods. 


where [_-Jn denotes the ra-th digit truncation. 

So we first construct a classical (order 1) digital net from the Xing-Niederreiter sequence 


using the ’sequence-to-net’ propagation rule (6.7). Then we employ the digit interlacing oper¬ 


ation (4.1) to obtain an order 2 net. 


For this particular kernel and point set the formula for the squared worst-case error can be 
written as 


sup 

ll/l^(nix( Td )H< 1 


Ifd 


2’ 1 —1 

f(x) dx - 2~ n /(y 

i =o 


2 n — 1 2 n —1 


_1 + 2 -2n XX K d Jx\: 


i =0 j =0 


( 6 . 8 ) 


In Figures [2] and |3] we computed the worst-case errors of the described construction of an 
order 2 digital net in H } n - ^(T d ) for dimensions d = 2,..., 5 and compared it to the bounds 
These expected rates of convergence A r_r (logIV)^ _1 ' >//2 were plotted in 


5.3 


from Theorem 

dashed lines. One can see that our observed rate of convergence matches the predicted one 
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even though a dimension-dependent constant seems to be involved. Additionally we computed 
the worst-case error for the Halton construction [23] which is amongst the most popular QMC 
sequences and performs very well for smoothness r = 1. Moreover, we consider the sparse 
grid construction which consists of certain tensor products of the univariate trapezoidal rule, 
yielding an error decay of 0(N~ r (log iV)C rf_1 )C T ’+ 2 )), see m • Sparse grids go back to ideas from 
Smolyak [35] and belong to todays standard approaches when it comes to high-dimensional 
problems, see e.g. |S] and the references therein. It can be seen that their rate of convergence 
depends stronger on the dimension d than the low-discrepancy approaches. 

The same analysis was done for the case of second order smoothness in H^ lix (T d ). The 
results are given in Figures [4] and [5] Here, all quantities were computed in 128-bit floating 
point arithmetic. In the bivariate case it is known [ 351 [20] that the Fibonacci lattice performs 
asymptotically optimal. This can also be observed in Figures [2] and [4] where the Fibonacci 
lattice yields the same (optimal) rate of convergence as the order 2 digital nets, although it 
seems to have a significantly smaller constant. For small Fibonacci numbers, it is even known 
that the Fibonacci lattice is the globally optimal point set [26] • In summary, we can see that 
the order 2 net, the Fibonacci lattice and the sparse grid are able to benefit from the higher 
order smoothness, while the Halton sequence does not improve over A' _1 (logIV)( rf ~ 1 ' ) / 2 . 


6.2 Integration of kink functions 


Mixed Sobolev regularity H^ ix is often not suitable to reflect the correct asymptotical behavior 
of the integration error of one fixed function. In case of kink functions, like for instance the 
Faber hat functions Vjj? from (3.2), we observe the Sobolev regularity Vjj~ G iL 3 / 2 ~ e whereas 
the Besov regularity is Bf^. The tensor product kink functions belong to H^- x £ , but as 
well to Sf^B. This can be easily deduced from the characterization in Lemma 2.7 Glancing 


at Theorem 5.3, we see that the (optimal) error bound does not depend on the integrability 
parameter p of the mixed Besov space S r p q B. Hence, it seems to be reasonable to “sacrifice” 


integrability in order to gain smoothness which makes our Besov model more suitable for this 
issue. Our first example is a typical kink function of the form g(x) = max{0, h(x)}. To be 
more precise, we consider tensor products of the univariate ( normalized) function 


g (x) = max | ^ - (x - 1/2) 2 ,0 


(6.9) 


which belongs to ^^(T) and has integral /[' g(x) dx = 1. Hence the tensor product function 

d 

9d(x) :=Y[9(xj) , xeT d , (6.10) 

j =l 

belongs to Sf 00 B( T cZ ) with integral f Jd gd(x) dx = 1 and the same holds for the shifted functions 

d 

9d(x , g) := JI g(frac(xj + rjj)) , x G T d , (6.11) 

j =i 

where frac(t) = t — [ij denotes the fractional part of t and g G [0, l\ d . 

In order to obtain smooth convergence rates we compute the maximum error of 1000 
randomly shifted instances of gd, be. 


Rn(9) — max 

fce{i,...,iooo} 


N 


/ T d 


9d(x , g k ) dx - ^ A mix 1 , g k ) 


i= 1 


( 6 . 12 ) 
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Figure 6: Maximum over the relative integration errors for the test function (6.11) with kinks. 


Here, the shifts i] k ~ U[ 0, l] d are independent and identically uniformly distributed in for 
k = 1,..., 1000 and the integration nodes x l ,i = 1,... ,N and associated integration weights 
A i depend on the chosen integration method and also the total number of function values N. 
The results are given in Figure [6j where we compared the performance of the order 2 nets to 
both the sparse grid and Halton construction. 

o /o 

Next, we consider a toy example from H 1 / C ^ ) (T) which has Sobolev regularity below r = 1. 
We take the square root of the level 0 hat function ( |3.2[ ) normalized with respect to Li(T d ), 
i-e., 

9 (*) : = ^\J v o,o{t) • (6.13) 


Hence the tensor product function 


It holds Jj g(t,) dt = 1. The Besov regularity r = 3/2 can be easily deduced from Lemma 


2.7 


d 

9d{x ) := Ylgixj) 

3 = 1 


x e T l , 


(6.14) 
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Figure 7: Maximum over the relative integration errors for the test function (6.15) with frac¬ 


tional smoothness r = 


belongs to Si^B(T d ) with integral f Td g d (x) dx 


1. The same holds for the shifted functions 


d 

g d (x, rf) := g(hac(xj + rjj)) , x <E T l . 
i=l 


(6.15) 


Again, we compute the maximum integration error of 1000 shifted instances of g. The results 


are given in Figure [7J It can be c 
the ones predicted in Theorem 


5.3 


early observed that the obtained convergence rates match 
i.e. N~i (log A r ) d_1 . 
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